O 
O 
(N 

Cu- 
te 

o 

> 

in 

m 
m 

o 

O 

Or 

6 



13 



Cosmological parameters from combining the Lyman-a forest with CMB, galaxy 

clustering and SN constraints 

Uros Seljak, 1,2 Anze Slosar, 3 and Patrick McDonald 4 

'Department of Physics, Princeton University, Princeton NJ 08544, U.S.A. 

2 International Center for Theoretical Physics, Trieste, Italy 

3 Faculty of Mathematics and Physics, University of Ljubljana, Slovenia 

^Canadian Institute for Theoretical Astrophysics, University of Toronto, ON M5S 3H8, Canada 

(Dated: February 3, 2008) 

We combine the Ly-a forest power spectrum (LYA) from the Sloan Digital Sky Survey (SDSS) 
and high resolution spectra with cosmic microwave background (CMB) including 3-year WMAP, 
and supernovae (SN) and galaxy clustering constraints to derive new constraints on cosmological 
parameters. The existing LYA power spectrum analysis is supplemented by constraints on the mean 
flux decrement derived using a principle component analysis for quasar continua, which improves 
the LYA constraints on the linear power. We find some tension between the WMAP3 and LYA 
power spectrum amplitudes, at the ~ 2<r level, which is partially alleviated by the inclusion of other 
observations: we find as — 0.85 ± 0.02 compared to a% = 0.80 ± 0.03 without LYA. For the slope we 
find n s — 0.965 ± O.Of 2. We find no evidence for the running of the spectral index in the combined 
analysis, dn/dlnk = — (f.5 ± f.2) x 10 -2 , in agreement with inflation. The limits on the sum of 
neutrino masses are significantly improved: Yl m v < 0.17eV at 95% (< 0.32eV at 99.9%). This 
result, when combined with atmospheric and solar neutrino mixing constraints, requires that the 
neutrino masses cannot be degenerate, mz/mi > f.3 (95% c.L). Assuming a thermalized fourth 
neutrino we find m s < 0.26eV at 95% c.L and such neutrino cannot be an explanation for the 
LSND results. In the limits of massless neutrinos we obtain the effective number of neutrinos 



m 



5.3 



+0.4+2.1+3.8 
-0.6-1.7-2.5 



and N£ = 3.04 is allowed only at 2.4 sigma. The constraint on the dark 



energy equation of state is w — —1.04 ± 0.06. The constraint on curvature is f2k = —0.003 ± 0.006. 
Cosmic strings limits are Gfj, < 2.3 x 10 -7 at 95% c.L and correlated isocurvature models are also 
tightly constrained. 



PACS numbers: 98.80.Jk, 98.80.Cq 



I. INTRODUCTION 

The latest WMAP results of cosmic microwave 
background (CMB) temperature and polarization 
anisotropies from 3 years of observations have presented 
a remarkably clear view of the universe in its early epoch 
around redshift 1100 0,|2|- These data, combined with 
other tracers such as galaxy clustering and supernovae, 
can be satisfied within a simple cosmological model in 
which the universe is filled with cold dark matter and 
dark energy and the initial conditions are generated 
through a process like inflation [3j . One of the main new 
results from this analysis relative to the 1st year WMAP 
analysis jj| is the reduced value of the optical depth due 
to reionization, which follows from a more sophisticated 
modeling and better statistics of CMB polarization data. 
Many other parameters are affected by this change, 
such as a reduced value of the amplitude of fluctuations 
a& and increased evidence for a red spectral index of 
primordial fluctuations, n s < 1 (note, however, that the 
optical depth constraint accounts for only part of the 
change in parameters). 

To increase the statistical power of constraints the 
WMAP team added other tracers, such as small scale 
CMB pjiBH, galaxyclustering from 2dF and SDSS 
yj [T(J, El i supernovae [Tj, [pj and weak lensing [lj, [l5| . 
All of the tracers used in their analysis trace large scale 



structure on scales above ten megaparsecs. However, 
there are many physically motivated models that can 
only be distinguished using small scale information. One 
is the running of the spectral index, a = dn/dlnk. Infla- 
tion generically predicts that there should be no running 
at the observable level in the current data. Therefore, 
any detection of running would be a major surprise and 
would force us to either accept models of inflation with 
an unusual shape of the potential with large third deriva- 
tive or abandon it in favor of a better theory. In the lat- 
est WMAP analysis the data are marginally better fitted 
with non-zero running, a = —0.06 ± 0.03, but this anal- 
ysis did not include small scale information, so the error 
is large. 

Another class of models where small scale clustering in- 
formation helps are those with massive neutrinos. These 
suppress power in dark matter clustering on small scales 
because of their free streaming, which erases their own 
fluctuations on scales below the free streaming length. 
This in turn slows down the growth of cold dark matter 
structure on the same scales, leaving an imprint in the 
dark matter power spectrum. For the relevant range of 
neutrino masses below leV the full extent of the suppres- 
sion happens on megaparsec scales. From atmospheric 
and solar neutrino mixing detections we know that at 
least one neutrino has mass above 0.05eV and another 
above O.OIeV [Ig, H3, but these experiments only tell 
us the mass difference squared and cannot provide the 



absolute scale for neutrino mass. Cosmological observa- 
tions weigh neutrinos and can provide an absolute scale 
for their mass. Thus cosmological observations can an- 
swer one of the most important issues in neutrino physics 
today, the nature of the mass hierarchy of neutrinos. 

It is clear that to improve upon these constraints one 
should determine the fluctuation amplitude on smaller 
scales than probed by CMB, galaxies or weak lensing. 
Nonlinear evolution prevents one from obtaining useful 
information at z = 0, so one must look for probes at 
higher redshift. Of the current cosmological probes, the 
Ly-a forest (LYA) - the absorption observed in quasar 
spectra by neutral hydrogen in the intergalactic medium 
(hereafter IGM) - has the potential to give the most pre- 
cise information on small scales 18]. It probes fluctua- 
tions down to megaparsec scales at redshifts 2-4, so non- 
linear evolution, while not negligible, has not erased the 
primordial information. 

Currently the most precise measurement of the LYA 
power spectrum comes from the analysis of more than 
3000 spectra from the Sloan Digital Sky Survey [iHI^. 
This data set is almost two orders of magnitude larger 
than available previously and has dramatically reduced 
the errors on the amplitude and slope of the power spec- 
trum around wave- vector k ~ 1/iMpc -1 . Combining 
this and a small amount of higher resolution Ly-a for- 
est data [HI |21| with other probes has already been 
demonstrated to significantly improve the constraints on 
the running and neutrino masses [22j. For the running 
the current constraint from WMAP 1st year and LYA 
is a = —0.003 ± 0.010, which is nearly a factor of 3 
improvement over the WMAP 3 year analysis without 
LYA. For neutrino mass the constraints including LYA 
are X)m„ < 0-3 - 0.4eV at 95% confidence level jlll^ 
in the most restrictive 6 parameter analysis. 

The previous joint LYA and WMAP first year analy- 
sis in 22] suffered a significant increase in the errors on 
the amplitude and slope of the power spectrum from the 
Ly-a forest due to projecting out all of the poorly known 
nuisance parameters such as the temperature-density re- 
lation, flux of UV radiation and its fluctuations, filtering 
length, effects of galactic winds etc. The largest part of 
this increase is due to a partial degeneracy between the 
mean flux decrement determined by the intensity of the 
UV background and the amplitude of the power spectrum 
|18L [l|| • While this degeneracy is partially broken by the 
flux power spectrum information itself, the error on the 
amplitude increases significantly over the case where the 
mean flux decrement was known. 

One can improve upon this by providing external con- 
straints on the redshift evolution of the mean flux decre- 
ment. In principle, if one knew the shape of the quasar 
spectrum, one would extract the decrement simply by 
comparing the actual flux level to the predicted flux level 
in the absence of absorption. One way to determine the 
continuum level of quasars is to identify regions with no 
absorption, but this is difficult at redshifts above 2.5, 
where the amount of absorption is significant and such 



regions are rare. Another way is to assume quasar spec- 
tra all have the same shape and the continuum is known 
from observations at lower redshifts where absorption is 
negligible. In practice the quasars do not all have the 
same shape, so this approach has to be generalized. One 
can use the fact that the quasar continuum follows well 
defined variations and, identifying these at rest wave- 
lengths longer than 1216A, where there is no absorption, 
predict the shape in the Ly-a forest. A formalization of 
this procedure, including the errors, can be done by ap- 
plying the principal component analysis (PCA) on the 
measured spectra, possibly including those measured at 
low redshifts where there is little Ly-a forest observed 
|24j . This procedure leads to a better determination of 
the mean flux decrement as a function of redshift up to an 
overall amplitude, which cannot be determined with this 
method. However, at low redshifts the level of absorption 
is sufficiently low that the continuum can be determined 
simply by connecting regions of no absorption, especially 
when using high resolution spectra |2l|. This method 
therefore provides external information on the mean flux 
decrement, which can be used in the joint analysis. 

The procedure described above has been applied to 
SDSS data and results presented in |ig. They show a 
remarkable agreement between the PCA and power spec- 
trum based determination of the mean flux decrement as 
a function of redshift. This agreement provides addi- 
tional confidence in the results obtained from the power 
spectrum analysis. Since the PCA method gives some- 
what smaller errors on the mean flux decrement one can 
use these as external constraints to better determine the 
amplitude of the power spectrum. In the current analysis 
the improvement is modest, 20-30% on both amplitude 
and slope, because we do not have a tight constraint on 
the overall normalization of the mean absorption. The 
purpose of this paper is to apply these LYA constraints 
to the joint cosmological analysis with other data sets, 
including the new WMAP 3 year data. An independent 
analysis of similar type has recently been performed by 



II. ANALYSIS 

We combine the constraints from the Ly-a forest 
(LYA) [ji| with the SDSS galaxy clustering analysis 
(SDSSgal) dO], SDSS luminous r ed g alaxy constraints 
on the acoustic peak (SDSS BAO) [llf, 2dF galaxy clus- 
tering analysis (2dF) 9], Gold sample of supernovae from 
[121 ] and SNLS supernovae sample [l3J , CMB power spec- 
trum temperature and polarization observations from 
WMAP 3 year analysis (WMAP3) 0,0 and from smaller 
scale experiments, including Boomerang-2k2 8], CBI 5], 
VSA @, and ACBAR (CMBsmall). Contrary to the 
WMAP approach we do not include the SDSS bias con- 
straints [26} nor weak lensing constraints [lj], because 
their statistical power is weak when the LYA power spec- 
trum is included and we prefer not to cloud the interpre- 



tation of the results. We mention however that both of 
these constraints prefer higher normalization of erg than 
the recent WMAP3 results and provide additional sup- 
port for the results found in this paper. Moreover, for 
the case of weak lensing, the recent detection of intrinsic 
alignment J27j may push the weak lensing normalization 
even further up from existing numbers, which do not ac- 
count for this effect. 

Our most general cosmological parameter space is 

p = (wb,wdm,fik,fl, t, /„, 7V° ff ,w, A, n s ,a,r, 

A iso ,GiJ,), (1) 

where Wb = Ob ft 2 , with Ob baryon density in units of 
the critical density, Wdm = Od m ft 2 with Odm dark matter 
density in units of the critical density, h is the Hubble's 
constant in units of 100 km s — 1 Mpc -1 , £\ = 1 — O tot , 
where O to t is the total density of the universe in units of 
the critical density, 9 is the angular scale of the acous- 
tic horizon, r is the optical depth to the surface of last 
scattering, f v is the fraction of dark matter in neutri- 
nos (the rest is assumed to be cold dark matter), N° s 
is the the energy density in relativistic background ex- 
pressed in units of massless neutrino families and w the 
equation of state of dark energy, (which in general can 
be time dependent, but in this paper we only explore 
it as a constant). Parameters A- lso and Gfi describe the 
iso-curvature modes and string contribution (see section 
IIII I|) . The power spectrum of primordial fluctuations is 
assumed to be of the form 



Pi(k) = A 



feo 



ris+a log(fc/fc )/2 



with the pivot point chosen to be fco = 0.05/Mpc. Wc 
use two additional amplitude parameterizations, which 
are derived parameters: a% is the rms of fluctuations in 
the linear density field smoothed at 8/i _1 Mpc scale at 
z = 0, and A 2 (fc = 0.009s/km, z = 3) describing the 
fluctuations at roughly one megaparsec scale at z = 3. 
Finally, r is the ratio of tensor to scalar power spectrum 
amplitude at pivot point fco. 

We begin the exploration with the simplest 6 parame- 
ters (p = LOh, Wdmj 6, T , A, n s ) and then proceed by adding 
additional parameters, one at a time. We have tested 
the analysis using two independent Monte Carlo Markov 
Chain (MCMC) codes. The first code is described in 
the previous analysis of the Ly-a forest J22| and uses 
CMBFAST 0]. The second is the standard CosmoMC 
package J2£| using CAMB [20]. We find a very good 
agreement between the two codes. We also compare 
our results whenever possible to the WMAP analysis, 
finding very good agreement in the analysis of WMAP 
data. We find some small discrepancies when other data 
sets are combined with WMAP. For example, in our 
WMAP+SDSSgal analysis we find ct 8 = 0.803 compared 
to 0.772 in the WMAP paper :3], which we suspect can be 
explained by the different treatment of nonlinear effects 
in the SDSSgal data, as discussed below. 



For galaxy surveys we adopt a somewhat more conser- 
vative approach in the modeling of nonlinear effects than 
what was done previously. The nonlinear biasing effects 
are argued to be well described by the expression given 
in |3, 



.Pgal(fc) = b' 



1 + Qfc 2 
1 + A g k' 



Pr, 



xW. 



(2) 



We use this expression with A g = 1.4 and marginalize 
over the bias parameter b, but rather than adopting a 
fixed value for Q we treat it as a free parameter. For the 
2dF power spectrum we adopt a Gaussian prior Q = 4.6± 
1.5, as found in |j|. For the SDSS galaxy power spectrum 
we use the same approach with Q = 10 ± 5 and Q > 0. 
We use a higher mean value for Q because the SDSS 
data were analyzed in real space, where nonlinear effects 
are stronger and the value Q — 10 fits well the standard 
nonlinear power spectrum from N-body simulations (at 
a percent level). Exploring the allowed range of Q in 
numerical simulations we found that even values above 20 
are possible in real space. Therefore, we argue that using 
a fixed value for Q is not justified given the uncertainties 
in the nonlinear galaxy biasing. Instead we allow for a 
wider range of Q to account for the uncertainty and let 
the data choose the best value. In both cases we use the 
data up to fc = 0.15/i/Mpc, so the nonlinear effects are 
small, but not negligible. 

For baryonic oscillations (BAO) we use the recent mea- 
surements of the baryonic peak in the Luminous Red 
Galaxies sample from the SDSS as published in [llj ■ To 
compare theoretical power spectra with the measured 
correlation function, we follow the prescription of the 
authors of the data: first we interpolate between the lin- 
ear power spectrum and the "no-wiggle" equivalent using 
the weighting function exp(— (afc) 2 ) with a = 7/i~ 1 Mpc. 
For the "no-wiggle" spectrum we took either the ap- 
proximation from [31| (as used in |ll|') or the convolu- 
tion with suitable top-hat in the logP(fc)-logfc space. 
This is because the approximation of |3JJ is not avail- 
able for all cosmological models. We also compared the 
results to those found when this step is ignored. In gen- 
eral we find that the effects from these different treat- 
ments on the final results are below one sigma, but when 
this is not the case we explicitly discuss it. Then the 
non-linear corrections are applied using the HaloFit [321 ] 
package. The resulting power spectrum is converted to 
the corresponding correlation function. Finally, the red- 
shift space distortion and non-linear bias effects are ac- 
counted for by multiplying the correlation function £(r) 
by (1 + 0.06/(1 + (0.06/i/Mpcxr) 6 )) 2 , in agreement with 
TV-body simulations. This correlation function is then 
compared with the measurement using the full covari- 
ance between errors. The amplitude of the correlation 
function is marginalized over. 
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B. Amplitude of fluctuations: 
compatible? 



are the data sets 



Wdir 

h 



0.117 



+o. 



003+0.005+0.007 
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0.7051^ 
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1.012+0.025+0.037 
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Ts 
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< 0.22 (< 0.37) 

n nl r+0.012 +0.023 +0.036 
U.U10_ 0012 _o. 021 -0.029 



n k 

w 
4- , 

^Mso.bar 

Ai so ,CDM 

Gfl 



5.3 



-0.4+2.1+3.8 



-0.6-1.7-2.5 

<0.17eV(<0.32eV) 
< 0.26eV (< 0.43eV) 



-0.003 



+0.0060 +0.0109 +0.0157 
-0.0061-0.0122-0.0180 



1 n/in+ ' 063 +°' 124 + ' 178 

l.Uy:U_ .063-0 .130-0 .208 

_n nfi+ 018 +°- 35 + - 50 

u - UD -0. 18-0. 34-0. 55 

_n nn7+°- 034 + - 068 + ' 110 

U.UUr_ . 035 -0.067-0.104 

< 2.3 x 10~ 7 (< 2.9 x 10~ 7 ) 



TABLE I: Parameter constraints for the simplest 6-parameter 
model. Parameters below the horizontal line (r. . . Gfi) were 
constrained individually, varying the basic 6 parameters and 
one extra parameter at a time. Limits are 1, 2 and 3 sigma; 
when only upper limits are given, these are at 95% and 99.9%. 



III. RESULTS 



A. Basic cosmological parameters 



Table Q] shows the results of the analysis with 
and without LYA for various parameter combinations, 
starting with the simplest 6 parameter model p — 
(t, Wb,Wdmi ^a, A, n s ). For most parameters adding LYA 
reinforces the conclusions derived by WMAP. However, 
LYA prefers a high value of the normalization and as a 
consequence other parameters arc also affected. For ex- 
ample, the preferred optical depth in the joint analysis 
increases to r = 0.11 and the normalization increases to 
er 8 = 0.85. 



The parameter where there is the most tension be- 
tween the analyses with and without the Ly-a forest is 
the amplitude of fluctuations. In terms of a%, in [22| 
this parameter was found to be <r$ — 0.90 ± 0.03, but 
this value was driven by the WMAP optical depth prior, 
as seen in Figure 1 of that paper, which suggests that 
for the new value r = 0.09 one would find as = 0.84. 
In our new analysis we find a% — 0.85 ± 0.02, which is 
in a good agreement with this expectation. This should 
be compared to cr§ = 0.75 ± 0.06 for WMAP alone and 
a s = 0.80 ± 0.03 for WMAP combined with everything 
else in our analysis. There is some tension between the 
WMAP and LYA values and the most relevant question is 
whether this can be explained as a statistical fluctuation 
or as a signal of a systematic problem in the data. 

The Ly-a forest effectively constrains only the ampli- 
tude and slope of the linear power spectrum at pivot 
point k = 0.009 s km~ and z = 3. Projections to 
other scales or redshifts are sensitive to cosmological pa- 
rameters that are essentially unconstrained by the Ly-a 
forest. On the other hand, WMAP alone can make a 
well-controlled projection into the LYA amplitude-slope 
plane, so we make our comparison in this plane. Con- 
straints from various experiments are shown in Figure ^ 
The top panel is for the case of no running using the ba- 
sic 6 parameters. In this case WMAP data alone already 
provide a strong constraint in this plane. We see that 
LYA data prefer a higher value for the amplitude than 
WMAP and about the same value for the slope, but the 
2-sigma contours have a considerable overlap. The agree- 
ment is improved if other probes, such as SDSS galaxy 
clustering, are added to WMAP data, which is also shown 
in the figure. In terms of differences in x 2 between fit- 
ting for CMB+SDSSgal+SN and LYA separately versus 
jointly we find that a joint fit increases the total \ 2 by 3.8, 
also suggesting less than 2-sigma discrepancy. Finally, we 
see that the Ly-a forest constraints using only high res- 
olution spectra |33l l34| , while considerably weaker com- 
pared to our analysis, also prefer a higher normalization. 
Here we have taken the values given in 35], including 
the 15% error on the amplitude, marginalized separately 
in two redshift bins. We note that our own analysis of 
these data gives a somewhat smaller error on the ampli- 
tude and a somewhat larger error on the slope, but the 
basic conclusions remain the same. 

The situation is similar when running is added as a 
parameter. One can see that WMAP3+SDSSgal+SN er- 
rors are comparable to LYA. In this case the WMAP3 
data prefer large negative running, so even though their 
errors significantly increase in this plane, the mean value 
has been pulled towards lower values for both amplitude 
and slope and the agreement with LYA is actually slightly 
worse than before. Still, there is considerable overlap of 
2-sigma contours. The difference in y 2 between fitting 
WMAP3+SDSSgal+SN and LYA separately and jointly 
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FIG. 1: Constraints on effective slope n e ff and amplitude A 
of the linear power at k p — 0.009 s/km, z — 3. The top 
panel corresponds to models with power-law primordial fluc- 
tuations, while the bottom includes running of the spectral 
index. In all cases 68% and 95% contours are shown. Red, 
short-dashed contours are for WMAP3, black, long-dashed 
for WMAP3+SDSSgal+ SN, green, solid for LYA only, and 
yellow, solid contours with horizontal stripes are for all of 
these constraints combined. The blue, dotted contours show 
an alternative Ly-a forest analysis using only high resolution 
spectra [351 ]. Top panel also shows the WMAP3 constraints 
when N„ is allowed to be free (thin red short-dashed) versus 
being fixed at N° s = 3.04 (thick red short-dashed), as well as 
the central values of WMAP3 contours for different values of 
Nf from 1 to 10. 



is now 6. 

So is there any evidence for inconsistency between 
WMAP3 and LYA data? From the comparisons made 
above it would seem that in the simplest 6-parameter 
models the discrepancy is around 2-sigma. However, it is 
important to recognize that any attempt to answer this is 
an a-posteriori attempt and can only provide qualitative 
answers. For example, if we focus on a single parameter 
like amplitude, which is chosen a posteriori from N pa- 
rameters on the basis of the fact that it shows the most 
discrepancy, then we will obtain biased results: compar- 
ing two experiments in N parameters and selecting the 
most discrepant one on the basis of one out of p crite- 
ria (for example, if one can choose from larger than and 
smaller than then p=2) will on average give one parame- 
ter in the 1/pN probability corner. In our case LYA most 
naturally measures amplitude and slope at one mega- 
parsec, so N=2 and p=2 and the statistical significance 
of discrepancy is further diluted by a factor of 4. Our 
conclusion is thus that there is no compelling evidence 
that the differences between LYA and WMAP3 cannot 
be explained as a normal statistical fluctuation and we 
proceed under this assumption. We return to this in the 
discussions. 



Testing inflationary models: 
tensors 



spectral index and 



LYA prefers a higher amplitude at a megaparsec scale 
and this also pushes the primordial slope upwards. The 
effect is small, but relevant for the question of whether 
scale invariant spectrum with n s — 1 is statistically 
excluded. We find n s — 0.965 ± 0.012 versus n s = 
0.951 ± 0.013 without LYA, so about one sigma change 
in the value. A very similar value was obtained in [36}. 

Adding tensors to the analysis does not improve the 
fit. We find 



r < 0.22at95%c.l. (< 0.37 at 99.9%), 



(3) 



slightly stronger than r < 0.28 — 0.3 at 95% c.l. from 
WMAP3+SDSSgal or WMAP3+2dF analysis % 

In Figure |21 we show the constraints in (n s ,r) plane. 
This plane is of interest for inflation model builders, since 
many of the inflationary models can be classified accord- 
ing to their predictions in this plane. We see that models 
with n s > 1 are disfavored at r = 0, but there is a positive 
correlation between r and n s . 



D. Testing inflationary models: running of the 
spectral index 

The issue of the running of the primordial slope con- 
tinues to be controversial. The latest WMAP analy- 
sis finds the strongest evidence for running so far, par- 
ticularly when combined with other data sets such as 
small scale CMB or SDSSgal, e.g. a = -0.066taol for 
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FIG. 2: This figure shows constraints on the n s -r plane. Con- 
tours enclose 68% and 95% probability. 



WMAP+CBI+VSA. The dynamic range of scales is lim- 
ited by CMB and galaxy clustering, and the correspond- 
ing error is large. It was previously shown that the error 
is significantly reduced when SDSS LYA and WMAP first 
year data were combined, giving a = —0.003 ±0.010 |22| . 

Here we redo the analysis using the new WMAP3 data 
and our LYA analysis. Together with all the other data 

we find 



-0.015 



-0.012- 
-0.012- 



0.023- 
0.021- 



0.036 
0.029 



(4) 



The combined analysis therefore continues to show no 
evidence for running, in agreement with theoretical ex- 
pectations from inflation. The errors are improved rel- 
ative to the WMAP 3 year analysis without LYA, but 
have not improved relative to LYA + first year WMAP, 
which also did not show any evidence for running. At 
this point, a conservative conclusion is that the results of 
the WMAP team analysis without LYA are a statistical 
fluctuation which is eliminated by adding the LYA data, 
which arc much more constraining than any combination 
without LYA. We note that 2-sigma deviations from the 
true value on one of a dozen parameters are bound to 
happen in almost any analysis and the above discussion 
of dangers of a-posteriori assessment of statistical signif- 
icance also applies here. 



E. Number of relativistic neutrino families 

In the standard model the main relativistic compo- 
nents of energy density are photons and neutrinos. Neu- 
trinos decouple earlier than photons, which receive addi- 
tional entropy transfer from electron-positron annihila- 
tion, and as a result photons are hotter than neutrinos, 
T v = (4/ll) 1/3 T 7 . The neutrino contribution to the en- 
ergy density is usually described in terms of an effective 
number of neutrino families assuming the temperature 
relation above. The standard model predicts N° =3.04 
instead of 3 because some of the entropy transfer goes to 
neutrinos that are not completely decoupled by the time 
of electron-positron annihilation. More generally, we use 
N^ s as a parametrisation of relativistic energy density 
in the early universe, although it should be noted that 
effects of neutrinos may differ from other relativistic par- 
ticles (33 • If we assume N° s neutrino families then the 
analysis finds 



\reS _ r q+0.4+2.1+3.8 
ly v — °'°-0.6-1.7-2.5 



(5) 



We find that P(Nf) < 3.04 is 8.5 x 10~ 3 . If only in- 
teger values of the number of neutrinos are allowed, we 
find that 3:4:5 neutrinos are allowed with relative 
probabilities 1 :~ 6.4 :~ 9.7. How robust is the evi- 
dence against 3 neutrino families? Adding the HST con- 
straint 38] in the form of a Gaussian prior on Hubble's 
constant h = 0.72 ± 0.08 lowers Nf to 4.8±g;f ±i:|±|;? 
(but leaves the lower 3-sigma bound largely unchanged) . 
This is because there is a strong degeneracy between 
the number of relativistic degrees of freedom and the 
dark matter energy density, as any additional relativistic 
component can be compensated by additional dark mat- 
ter in order to keep matter-radiation equality unchanged 
[33 and so strong degeneracy exists with h and uJdm- 
For similar reasons, the constraint is also partially sensi- 
tive to the BAO treatment. Without BAO we find that 



N' 



6.0 



+1.4+2.9+3.8 

-1.4-2.4-3.3- 



Removing Lyman alpha data lowers the number to 
^■^-o.t-i'/r-i'Jn so tne evidence is significantly weakened. 
However, this is with both SDSS and 2dF galaxy power 
spectra and there is significant variation between the two. 
We find that the WMAP3 data together with super- 
novae data and SDSS main sample galaxy data favour 
high Nv = 7.8±J^±|;|±|:|, while replacing SDSS galaxy 
data with 2dF data favours much lower number of neutri- 
nos: Nf* = 3.2:^:0-2:3-3.0 (note that we cannot repro- 
duce the results in Table 10 of Q, however, the implau- 
sible asymmetry of their errors suggests some problem 
with those results). Excluding both the 2dF and SDSS 
main galaxy samples from the mix of all data results 
in 7V® ff = 5.2^06^1 8^2 6' wruc h is ornv a minor change 
compared to all datasets. 

We can conclude that the unexpectedly large result 
for the number of neutrino families does not rely on a 
single experiment as several combinations of data favour 
high values of relativistic energy density, albeit with an 
expected lower statistical significance. However, at the 



moment the evidence is strongest with the Lyman alpha 
data. Top panel of figure ^ shows why relaxing N° s im- 
proves the fit. When the parameter increases it leads to a 
higher amplitude and slope in LYA plane relative to the 
standard value, reducing the discrepancy between LYA 
and WMAP3 data. This is the only parameter we have 
found that has this effect. 

Our result of combining all datasets disfavors just three 
mass-less neutrinos, favouring higher numbers of neutri- 
nos instead. Big Bang Nucleosynthesis (BBN) in gen- 
eral tends to favour a lower number of neutrinos. In 
particular, ;4Q] obtain Nf = 3.14+g;™ from BBN data 
alone and Nf = 3.24+%%) by using the CMB photon-to- 
baryon ratio, while |41j get limits from N° s < 3.3 at 95% 
c.l. to A° ff < 6 and 95% c.l. depending on the details 
of the analysis. At the moment iV° ff = 4 is the inte- 
ger number of neutrinos most compatible with the data. 
Higher A° ff is at some tension with BBN constraints, 
while lower N® s has trouble with the datasets discussed 
in this paper. However, newly discovered systematic is- 
sues in various datasets might still shift values around. 
Finally, the possibility of a statistical fluke should not be 
discounted, as the standard value is within 2.4 sigma of 
the measured value. 
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FIG. 3: This figure shows the masses of individual neutri- 
nos plotted against the sum of masses which is constrained 
by the cosmological data. Solid line (red) is for normal hi- 
erarchy, while dashed (green) is for the inverted hierarchy 
(note that two upper to mass eigenstates are nearly degen- 
erate in this case). Vertical dotted lines correspond to our 
68%, 95% and 99.9% confidence limits. The assumed mass 



squared differences were taken to be AmJ 2 
and Am| 3 = 2.5 x l(T 3 eV 2 . 



x lO^eV 



F. Massive neutrinos 

As mentioned in the introduction, small scale cluster- 
ing combined with large scales and CMB is a very efficient 
way to weigh neutrinos, since these suppress the growth 
of structure in CDM. One can expect that adding LYA 
leads to a significant improvement in the constraints. 
Since the Ly-a forest prefers a high normalization com- 
pared to the CMB and the two data sets are barely com- 
patible with each other for massless neutrinos, adding 
massive neutrinos which suppresses power on the LYA 
scale would further increase the discrepancy and make 
them statistically incompatible, thus providing a tight 
limit on the neutrino mass. We find 






< 0.17eVat95%c.l. (< 0.32eVat99.9%). (6) 



These results depend somewhat on the treatment of the 
BAO data. If smoothing of the spectrum is used the 
limits raise to 0.35eV at 99.9% c.l. and to 0.40eV if 
BAO data are not included at all. Either way, this re- 
sult appears to rule out the claimed Heidelber g-M oscow 
experiment of neutrino- less double beta decay [42J ] . 

If we combine the constraint on the sum of the masses 
with the existing neutrino mixing constraints, which are 
Am 2 12 = 8 x lCT 5 eV 2 for solar and Am 2 3 = 2.5 x 10~ 3 eV 2 
for atmospheric mixing 16, 17] we find for a regular hi- 
erarchy mi ~ rrt2 < 0.05eV and 777,3 < 0.07eV, while 
for 2-sigma lower limit Ato 2 3 = 1.8 x 10 _3 eV 2 one finds 
to 3 < 0.065eV. At 2-sigma the degeneracy is thus lifted 
at the 30% level, 777,3/777,1 > 1.3. For an inverted hier- 
archy one has mi ~ rri2 < 0.062eV and 7773 < 0.045eV 



for Ato 2 3 = 1.8 x 10 _3 eV 2 and again the degeneracy be- 
tween mi and 771.3 is lifted at the 40% level with 2-sigma 
confidence. This is illustrated in Figure |3| 

It is interesting to note that with another factor of 
2 improvement we will be able to distinguish a normal 
hierarchy (for which ^2 fn v ~ 0.05eV) from an inverted 
hierarchy (with ^2 fn v ~ O.leV), even in the limit of one 
massless neutrino. This is also shown in Figure |3 where 
the 68% confidence limit is already below the minimal 
inverted hierarchy solution and an inverted hierarchy is 
disfavored at 80% confidence. We cannot yet differentiate 
between these two cases in the data with any statistical 
significance. We also note that assuming one or the other 
neutrino mass makes very little difference in the estimates 
of other data as there are no strong degeneracies with 
other parameters present. 

In addition to the strong evidence for mixing from 
solar and atmospheric neutrino experiments there is a 
hint for another oscillation with a much larger mass- 
squared difference coming from the LSND experiment 
|43j. This requires a mass squared difference of around 
leV 2 , which cannot be incorporated in the standard 3 
neutrino scheme. Of all the schemes that could incor- 
porate this result the one that is most plausible is 3+1, 
three active and one sterile neutrino. However, even this 
scenario is strongly constrained by the short baseline os- 
cillation experiments and the window is nearly closed at 
3-sigma level [13 • Existing cosmological constraints sug- 
gest that a sterile neutrino with mass at or above 0.8eV 
is ruled out at 2-sigma, but not yet at 3-sigma [3ij. So 
it is worth revisiting this issue with the new analysis. 
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Our new analysis considerably strengthen these limits. 
We find m„ < 0.26eV at 95% c.l. (0.43eV at 99.9%). 
This completely rules out a thermalized sterile neutrino 
with mass around leV as the explanation for LSND ex- 
periment. However, it is possible to suppress the produc- 
tion of sterile neutrinos in the early universe, for example 
by large lepton asymmetry (44j |. Therefore, there exists 
a possibility to accommodate these results with LSND 
by postulating that the neutrinos are not thermalized 
and their concentration is significantly suppressed rela- 
tive to a thermal distribution. Note that even though the 
LSND-kind of sterile neutrino is disfavoured, the data in 
general prefer more than 3 standard neutrinos worth of 
relativistic species in the early universe. Thus a light 
sterile neutrino may well be preferred by the data, as 
might be other light particles. 



G. Curvature 

Curvature is strongly constrained by CMB data be- 
cause it determines the angular position of the acous- 
tic feature through its effect on the angular diameter 
distance to the last scattering surface. However, other 
quantities such as the cosmological constant and mat- 
ter to radiation ratio also have a similar, albeit less pro- 
nounced effect, so the CMB alone cannot determine these 
parameters separately with a high accuracy. Adding ad- 
ditional information, such as the baryonic acoustic hori- 
zon, Hubble constant or redshift-distance relation from 
supernovae, improves the constraints considerably. In 
our analysis we find 



ft k = -0.003 



+0.0060+0.0109+0.0157 
-0.0061-0.0122-0.0180 



(7) 



when curvature is added to the basic 6 parameters. This 
is the strongest constraint on curvature to date and the 
data continue to show no evidence for it. It should be 
noted that the curvature constraints are affected by the 
choice of parameters and adding some, like a time de- 
pendent equation of state, can weaken these limits con- 
siderably. On the other hand, inflation as the leading 
paradigm for structure formation predicts no curvature, 
so there is good justification to ignore it in the remaining 
analyses. 



H. Dark energy 

Dark energy affects the rate of growth of structure, es- 
pecially for z < 1 where the dark energy is dynamically 
important. We have assumed in the analysis so far that 
the dark energy equation of state is w = — 1. Here we 
relax this and explore the constraints on w. In this paper 
we only explore the simplest case of a constant equation 
of state. It has been shown in [22J that the constraints 
assuming a constant w are very similar to constraints ob- 
tained for time dependent w in a 3-parameter expansion 



of the equation of state as a function of expansion fac- 
tor, assuming these are given at redshift 0.2-0.3, the pivot 
point for equation of state information. So our constraint 
on w is more general in the sense that it remains valid 
for smooth time dependencies of the equation of state, 
assuming it is interpreted as the value at z=0.2-0.3. We 
find 



w — _i n/in+ 0063 + 0124 + - 178 

w — 1 ' u ^ u -0. 063-0. 130-0. 208' 



(8) 



This is an improvement over other recent constraints and 
we see that w = — 1 remains very close to the best fit. 
Note that we have included perturbations in the dark 
energy both for w < — 1 and w > — 1 , although this makes 
very little difference in the constraints given how close 
the results are to w = — 1, where perturbations vanish. 
Most of the improvement comes from the reduced error 
in the amplitude of fluctuations provided by the Ly-a 
forest, which helps when combined with other data sets. 
Our results prefer somewhat more negative values for w 
than the WMAP analysis [3|, closer to the cosmological 
constant w = — 1. As a result those tracker quintessence 
models which predict w ~ —0.7 [45j appear to be strongly 
excluded. Other dark energy models which predict w ~ 
— 1 remain acceptable. 

It is interesting to investigate the sensitivity of the 
constraint to the inclusion or exclusion of supernovae 
data. Our results above include both supernovae from 
the SNLS sample [13| and from the gold sample of [l2fl . 
Without the latter we find w = —1.05 ±0.07, while with- 
out the former we find w = — 1.01 ± 0.07 and the two are 
statistically compatible. Finally, without any supernovae 
data at all the constraint is w = —1.02 ± 0.12, so the su- 
pernovae reduce the error by a factor of two relative to 
the other data. In all cases cosmological constant is the 
preferred solution lying within one statistical deviation 
of the best fit value. 



I. Iso-curvature and cosmic strings 

Finally, we constrain parameters of two possible exten- 
sions to the inflationary paradigm. Firstly we constrain 
the baryon or CDM iso-curvature modes, assuming com- 
plete correlation (or anti-correlation) and equal scale de- 
pendence of perturbations. Parameter A- lso is the ratio 
of amplitude in these fully correlated baryon/CDM iso- 
curvature modes to that in the adiabatic modes. We find 



-Ais. 



.4 



— _n nR+ 018 +°- 35 + - 50 

iso : bar — u - uu -0. 18-0. 34-0. 55 

n nn7+ - 034 +°- 068 + - 110 

iso,CDM - — U.UUf_ . 035-0. 067-0. 104- 



(9) 
(10) 



This is about a factor of 2 better than previous limits 
|46l |. For models predicting baryon isocurvature pertur- 
bations correlated with scalar perturbations from infla- 
tion the predictions for Aj SOj bar are to be of order unity 
[47| , which is severely constrained by the data. Curvaton 
model induced Aso,cdm & ve a larger range of possibili- 
ties, some of which remain viable even with the new data 
and corresponding analyses 46] . 



We also constrain the amount of cosmic strings allowed 
by the present data. We predict a cosmic string power 
spectrum using calculations in [4a, |43, Isjj , for a set of 
their fiducial cosmological and string network parame- 
ters. The only free parameter is the string tension G\x 
which accounts for a simple quadratic scaling in ampli- 
tude. We find that the present data require 

G^<2.3xl(T 7 at95%c.l. (< 2.9 x 1CT 7 at 99.9%). (11) 

This is comparable to other recent limits |49l l5lj . but 
is subject to several uncertainties in the string network 
evolution, which may affect the relation between string 
tension and CMB amplitude. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we combine the current Ly-a forest con- 
straints from SDSS, supplemented by a small amount of 
high resolution data, with the constraints from other cos- 
mological tracers, such as galaxy clustering, supernovae 
and CMB, including the recent 3 year WMAP analy- 
sis. Our Ly-a forest analysis is an improvement over the 
analysis presented in |19|. The main change is the in- 
clusion of the mean flux constraints (presented in |l9j 
but not used in the main power spectrum analysis) de- 
rived from a principal component analysis of quasar spec- 
tra, which allows one to identify the continuum level of 
quasars up to an overall amplitude by using the informa- 
tion from about 10,000 quasar spectra. These constraints 
are in remarkable agreement with those derived from the 
power spectrum analysis, providing additional support 
for the consistency of the overall picture of the Ly-a for- 
est. When the two are combined the resulting constraints 
on the amplitude and slope improve by about 20-30%. 

When adding this Ly-a forest information to other 
probes we find significant improvements on many of the 
parameters. For the running of the spectral index we 
find a = -0.015 ± 0.012, compared to a = -0.06 ± 0.03 
without the Ly-a forest. Our new value is in a good 
agreement with the previous analysis prior to WMAP3 
data [23] and in both cases provides no support for run- 
ning of the spectral index being non-zero. Our results 
are also in agreement with a similar analysis in J25J . 

The most dramatic improvements are found for neu- 
trino masses. This is expected, since the Ly-a forest is 
one of the few tracers of small scale power that is able 
to determine the amplitude of fluctuations in dark mat- 
ter on small scales, where massive neutrinos suppress the 
growth of structure. For the first time the constraints are 
sufficiently strong that, when combined with atmospheric 
and solar neutrinos, they lift the degeneracy of neutrino 
masses at the 30% level. The current constraints are 
beginning to distinguish between a normal and inverted 
hierarchy in the limit of zero mass for the lowest mass 
neutrino, with the inverted hierarchy being disfavored at 
80% confidence. Finally, we see no evidence for a fourth 



massive neutrino that would explain the LSND exper- 
iment. The upper limit on its mass is 0.26eV at 95% 
confidence limit, which rules out the entire range of pa- 
rameter space for the LSND experiment, assuming these 
neutrinos were thermalizcd in the early universe. On the 
other hand, a massless thermalized sterile neutrino is fa- 
vored, since we find Nf = 5.3±g;|±f ;^±1;|. 

Adding Ly-a forest data to the mix favors the sim- 
plest possible model without running, massive neutrinos, 
strings, isocurvature perturbations, curvature or dark en- 
ergy equation of state different from -1. The only ex- 
ception is the number of relativistic degrees of freedom, 
where the data prefer additional contribution beyond 
photons and 3 neutrino families. In many cases LYA 
brings the value of the parameters closer to the expected 
value than the analysis without LYA. As discussed in this 
paper, there is some tension between the Ly-a forest and 
CMB in the amplitude of fluctuations, but the statistical 
significance is small and the two data sets are compatible 
in amplitude at ~ 2-sigma. Any statistical assessment of 
compatibility must also take into account the a-posteriori 
nature of the comparison, since the data are measuring 
several parameters and one should not focus only on the 
most discrepant one to assess compatibility. It is worth 
mentioning that the amplitude tension of the two data 
sets is at a similar level of significance as the evidence for 
running or curvature in the WMAP data, which many 
(especially inflation theorists) would be willing to accept 
as a statistical fluctuation until proven otherwise. This 
is expected: one is estimating a dozen parameters and it 
is not surprising if some are 2-sigma away from their true 
values. In fact, the evidence for running is eliminated af- 
ter Ly-a forest data is added and the error is reduced by 
a factor of 2-3. So one can either accept that the original 
evidence was a statistical fluctuation which was settled 
conclusively with Ly-a forest and other data, or one can 
argue that some of these data have issues with system- 
atic errors and should be excluded from the analysis. In 
this paper we have proceeded under the assumption that 
these are just normal statistical fluctuations and that the 
data are compatible with each other. Future observations 
and investigations of possible systematic effects will pro- 
vide a more definitive answer. 

We wish to conclude with a brief discussion of the role 
of unaccounted systematics in the current cosmological 
analyses. No cosmological probe is completely immune 
to systematic uncertainties and all must be subject to 
ongoing investigations, as well as cross-checks among dif- 
ferent data sets. The remarkable progress in statistical 
precision of cosmological data over the past few years has 
been unprecedented, but often the improved statistics of 
the data sets has not been adequately matched by cor- 
respondingly improved analysis of systematic errors. To 
some extent this is expected: systematic errors are more 
difficult to track down and quantify and a detailed anal- 
ysis of these can take years, so is often done after the 
initial results have been published. However, such analy- 
sis is necessary and provides important confidence in the 
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results. A lot of progress has been made on this front 
in the past years and in some cases these concerns have 
been lessened after a detailed study of systematic errors 
has been published. A case in point is the remarkably 
detailed analysis of systematics in the recent WMAP 3 
year data, which has strengthened the reliability of cos- 
mological constraints derived from WMAP data. It is 
essential that corresponding investigations are done with 
other tracers as well and we believe that there is some 
room for improvement in many cases. For some cosmo- 
logical tracers, such as supernovae and weak lensing, the 
open issues (e.g. intrinsic alignments and shear calibra- 
tion in weak lensing, calibration and evolution in SN) are 
known and well organized campaigns to address them are 
underway. For other tracers, such as galaxy clustering, 
uncertainties in galaxy formation already limit the sta- 
tistical power that can be extracted on small scales, but 
the full extent of this uncertainty needs to be quantified 
better. 

Since the key results of this paper depend crucially on 
the reliability of the Ly-a forest it is particularly impor- 
tant to address the systematics of this tracer. Recent 
investigations of the Ly-a forest have shown that possi- 
ble astrophysical complications are either small or can be 
incorporated into the analysis J52, HE LlI uM ■ Moreover, 
an independent analysis of the SDSS flux power spectrum 
obtained results that are in agreement with the original 
analysis 1561 , although a somewhat lower amplitude is 
favored [25} (note that one does not expect exactly the 
same result because J5g used only SDSS data, while we 
add high resolution data that helps break parameter de- 
generacies). However, there are remaining worries about 
possible other effects that have not yet been thought of 
and more investigations are needed to exhaust the list of 
possibilities. 

One area where progress in Ly-a forest is expected 



soon is in numerical simulations. The original analysis in 
[jjj used an indirect approach patching together several 
simulations with differing box sizes to cover the whole 
dynamical range in observations and relying on hydro- 
PM simulations calibrated by only a few fully hydrody- 
namic simulations to cover the relevant range of param- 
eter space. With the fast hydrodynamic codes optimized 
for the Ly-a forest that are available today, such as Enzo 
|57j or the Eulerian moving frame TVD+PM code de- 
scribed in J53| , this is no longer necessary and parameter 
space can be covered entirely with hydrodynamic simu- 
lations. In addition, a code comparison project focusing 
on Ly-a forest is underway and will address the issue of 
agreement among different hydrodynamic codes. 

We expect that as the progress in statistical preci- 
sion of cosmological constraints slows down in the fu- 
ture more emphasis will be given to the investigation of 
systematic errors. Until then the current cosmological 
constraints obtained from combining many different cos- 
mological probes, including the ones presented in this pa- 
per, should be viewed as preliminary. However, the fact 
that the different tracers agree with each other within the 
range of normal statistical fluctuations bodes well for the 
current round of analyses and suggests systematic effects 
may be subdominant. 
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